Spin and charge dynamics of stripes in doped Mott insulators. 
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We study spin and charge dynamics of stripes in doped Mott insulators by considering a two- 
dimensional Hubbard model with N fermion flavors. For N = 2 we recover the normal one-band 
model while for N — » 00 a spin density wave mean-field solution. For all band fillings, lattice 
topologies and N = An the model may be solved by means of Monte Carlo methods without 
encountering the sign problem. At N = 4 and in the vicinity of the Mott insulator, the single 
particle density of states shows a gap. Within this gap and on rectangular topologies of sizes up to 
30 x f 2 we find gapless spin collective modes centered around q = (n ±e x ,n ±e y ) as well as charge 
modes centered around q = (±26^, ±2e y ), q = (ie x , ±e») and q = (0, 0). e WtV depends on the lattice 
topology and doping. 

PACS numbers: 71.27.+a, 71.10.-w, 71.10.Fd 



The understanding of the interplay between spin and 
charge degrees of freedom in doped two-dimensional Mott 
insulators remains a challenging issue. On the analyti- 
cal front, it is straightforward to apply the Hartree-Fock 
(HF) approximation but the difficulty lies in taking quan- 
tum fluctuations into account. On the other hand numer- 
ical methods take into account fluctuations but are lim- 
ited to small clusters due to the size of the Hilbert space, 
for exact diagonalization, or due to the minus sign prob- 
lem inherent to quantum Monte Carlo (QMC) methods. 
In this letter, we introduce a QMC method which lies 
between HF and full quantum fluctuations and apply it 
to the doped two dimensional Hubbard model. Since the 
TV— flavor model we consider breaks SU (N) spin symme- 
try it favors spin states. In particular in the vicinity of 
the Mott insulator (MI) we find a stripe phase experi- 
mentally observed in cuprates and nickelates jlj Our 
simulations reveal the dynamics of this phase. 

Our starting point is the Hamiltonian: 
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where i labels the sites of a square lattice and the 
first sum runs over nearest neighbors. The spinors 

el = ( ct • • • cl ] correspond to fermions with N fla- 

i \ i,l i,N J r 

vors. For even values of N, A Q . 7 = 5 a ^f(a) with 
f(a) = 1 for a < N/2 and -1 otherwise.' At TV = 2 
the model reduces to the standard 5T/(2)-spin invari- 
ant Hubbard model since the interaction is the square 
of the magnetization. Away from N = 2 the model has 
an SU(N/2) ® SU(N/2) symmetry which becomes clear 
when writing the Hamiltonian in terms of the spinors: 

(cl ■ ■ ■ ct , r ,„ ] and f ct r . ■ ■■ cl V As N -> 00 the 

V i,l i,N/2J V i, N/2+1 i,N ) 

exact solution of the model reduces to a spin density wave 
mean-field approximation. 



The central observation is that at N = An the minus 
sign problem in QMC approach is never present regard- 
less of the lattice topology and band-filling. To illustrate 
this and to keep the notation simple we consider the finite 
temperature auxiliary field QMC approach. After the 
usual decoupling of the interaction term with a Hubbard- 
Stratonovich transformation, the partition function be- 
comes Z = JV<j)e- NSW , where 
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H(r) 
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Here T corresponds to time ordering and in terms of 
our original fermions we can identify: cj + = c^ 1 and 
cv_ = c^jv/2+i' As — > 00 the integral over the 
fields (j) is dominated by the saddle point configuration, 
4>* satisfying dS{4>)/d<j)\^ =< i > ' — 0. A time independent 
solution is $* = (n? , — n^_) where the expectation 
value is taken with respect to the Hamiltonian in Eq. 
||. This is precisely the set of self-consistent equations 
obtained from the mean-field decoupling: — m t — — 
<f>* + [(rii. + — rii-) — $*] for the Hamiltonian of Eq. |^ at 
N = 2. 
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Tc 



is a real number for U > Q, it 

follows that NS(Q) is real for N — An. This observation 
allows us to interpret e - NS (®) as a probability distribu- 
tion which we sample with Monte Carlo methods. Hence 
for N = An the sign problem never occurs regardless of 
doping and lattice topology. For the simulations, we have 
used the related projector auxiliary field QMC algorithm 
which is better suited for the study of ground state prop- 
erties ||. Dynamical information is obtained by using 
the Maximum Entropy method 0]. The details of the 
approach will be discussed elsewhere. 
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FIG. 1: Ground state energy as well as spin and charge 
structure factors at q — (tv,tv). The data points at TV = 2 
stem from exact diagonalization studies Id]. 



To test and interpret our approach, we show in Fig. [I] 
the static spin and charge structure factors, 

S s (q) =^e i ^((no,+ ±n ,-)(n r , + ±n ri _)), (3) 

r 

as well as the ground state energy as a function of 1/N 
on a 4 x 4 lattice, U/t — 4 and two holes doped away 
from half-filling. The electron density is defined as: 
(n) = y ^2i( n i + + n i -} where V is the number unit cells 
i. In spite of sign problems for N ^ An, the QMC data 
interpolate between the mean-field solution at N = oo 
and the exact diagonalization results at N — 2. In the 
context of the large-N approach (|| Gaussian fluctuations 
correspond to 1/N corrections. Fig. [I] shows that this 
approximation is valid till N ~ 16. Hence, N — 4 is far 
beyond this approximation. In particular in the strong 
coupling limit, the model at N — 4 generates a mag- 
netic scale proportional to t 2 /U; a phenomena which is 
beyond Gaussian fluctuations around a weak coupling 
saddle point. Since we are solving the model exactly for 
a given value of N lattice symmetries, which are broken 
at the HF level, are restored in the QMC calculations. 

We use the method to investigate the metallic phase 
in the vicinity of the MI. To maximize quantum fluctua- 
tions, we concentrate on the N = 4 case. In the vicinity 
of the MI state, the relevant length scale is set by the 
inverse doping 1/8. To achieve this length scale at least 
along one lattice direction, we will consider rectangular 
topologies of width ranging up to 12 lattice constants. 
We adopt periodic boundary conditions in both lattice 
directions. 

We start with single particle excitations. The density 
of states N(oj) as a function of doping is plotted in Fig. 
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FIG. 2: Density of states as a function of doping. The lattice 
sizes are: (a)-(b) 28 x 8, (c) 30 x 8 (solid line) 30 x 12 dashed 
line, (d) Single particle occupation number. We have set 
U/t = 3, N = 4 and T = 0. 
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FIG. 3: Spin (a) and charge (b) structure factors as a func- 
tion of doping at U/t — 3. 



^ for various fillings. At U/t = 3, N — 4 and half- 
filling, we see a Mott gap of approximately 2t (Fig. ||a). 
Upon doping, the chemical potential shifts into the lower 
Hubbard band and spectral weight is transfered from the 
upper Hubbard band to generate low a energy feature 
0. In particular, at 8 = 1/14 remnants of the upper 
Hubbard band (« It away from fi) are seen and a low 
energy feature is detectable. The data of Fig. || also 
show the presence of a small quasiparticle gap at finite 
dopings. This is confirmed by the single particle occupa- 
tion number n(k) which shows a smooth behavior up to 
8 = 1/5. 

The origin of the quasiparticle gap lies in magnetic flue- 
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FIG. 4: Leading edge of the spin (a) and charge (b) dynamical 
structure factor as as function of width L y at U/t = 3, 8 = 
0.2, T = and TV = 4. We consider L x = 60 for L y = 1 
and = 30 otherwise, 

component of the wave vector has no meaning 



tuations which remain substantial away from half-filling. 
Fig. H shows the spin and charge structure factors for sev- 
eral dopings S on rectangles of width L y = 8. For those 
topologies, the dominant features are along the (1, 0) di- 
rection. In particular we see a peak in the spin structure 
factor at Q s = (n — e x ,7r) accompanied by one in the 
charge structure factor at Q c = (2e x ,0) with e x = ix8. 
Transforming the data in real space yields the stripe car- 
icature: rivers of holes along the y-direction separated 
by 1/8 lattice constants in the ^-direction. Across each 
river of holes, there is a 7r-phase shift in the spin structure 
with respect to the antiferromagnetic ordering. Similar 
features are seen in one dimension (ID) where our sim- 
ulations at N — 4 show dominant cusps in the charge 
(spin) structures factors at Akf = 2tt8 (2k / = tt — ttS) 
||. In ID we checked successfully for power law decay 
of the 2k f spin and 4fc/ charge correlation function at 
N = 4. On fixed width topologies one would equally ex- 
pect power law decay of the correlation functions. Wc 
were however unable to confirm this numerically due to 
limitations on lattice sizes. 

To study the dynamical properties of the striped phase, 
we compute the imaginary time displaced spin and charge 
correlation functions : S c i s (q, r) . By fitting the tail 
S c / s (q,T) to a single exponential, we obtain the leading 
edge of the charge and spin excitation spectra: w c / s (q). 
Irrespective of lattice width the overall features of the 
QMC data of Fig. || suggest gapless spin [charge] excita- 
tions with linear dispersions around Q s [Q c and (0,0)]. 
The data equally show other features, (i) At L y = 8, 12, 
the spin and charge excitations lie below the particle-hole 
continuum which from the single particle density of states 
(see Fig. ||(c)) lies at a threshold energy 2A qp « 0.5i. (ii) 
The continuity equation, 4zn(x, t) + V • J(x, t) = yields 



an identity between the dynamical charge structure fac- 
tor and the real part of the optical conductivity o~'(q, lv) : 
u>S c (q,u>)/q x = a' xx (q,u>). Thus long wave length charge 
modes carry current so that we have a metallic state. Di- 
rect QMC evaluation of the conductivity confirms this. 

It is tempting to try to understand the above exci- 
tation spectrum for our finite width lattices as evolving 
from ID physics ]8|, M. Here we would however like to dis- 
cuss two other possibilities in terms of Goldstone modes, 
(i) As mentioned previously, the model at N = 4 has 
an SU(2) ® SU(2) symmetry and one may be tempted 
to interpret the low lying spin and charge features in 
terms of Goldstone modes associated with the breaking 
of this continuous symmetry. Our calculations do not 
support this point view since those modes remain disor- 
dered, (ii) To interpret the data in terms of the dynamics 
of spin and charge density waves |l(], [ll[] , we have to ar- 
gue that the pinning of those modes due to the lattice 
is small. To do so, we carry out a HF calculation with 
order parameters: — n f.-) = A s cos(Q s f + <fi s ) and 

(rv,+ + n?,-) = cos(Q c r + <fi c ). At 8 = 0.2 we have to 
M Q s = G with G a reciprocal lattice vector and M = 10. 
The magnitude of the pinning decreases as a function of 
growing values of M and in the incommensurate case 
(M — ► oo) vanishes. To estimate the magnitude of the 
pinning, we solve the mean-field self-consistent equations 
for the moduli of the order parameters at fixed values of 
the phases. The mean-field energy then follows the form: 
E HF = E* F + A lC os 2 (0i + vr/2) + A 2 cos 2 ((/) 2 ) where 



= 20c 



and So = —c 



2<j) s . For the L y — 8 and 



L x = 30 system at 8 = 0.2 we find that A\ ~ — 5 x 10~ J 
while Ai ~ 0.8i. This implies that the pinning of the 6\ 
mode is extremely small since it is proportional to y ' A\. 
Following this point of view for the interpretation of the 
data, we will have to assume that the mass of the $i 
mode is beyond our numerical resolution so that we are 
not able to distinguish it from a genuine Goldstone mode. 
Pinning 62 = tt/2 and recalling that 6\ = 26 c + <f> s , and 
2Q C = Qs locks together the dynamics of the spin and 
charge density wave. The QMC results of Fig. | sup- 
port this point of view since the velocities of the spin 
and charge modes are comparable. 

We have up till now concentrated on narrow rectan- 
gular lattices where the topology forces stripes to occur 
along only one direction. As the width of the system 
increases to approach a square topology, peaks in the 
structure factors are seen along both lattice directions. 
Fig. [|a,b plots the spin and charge structure factors for 
our largest system: L y = 12, L x = 30. The spin shows 
pronounced features not only at (w — e x , w) but also at 
(71-, 7T — e y ) where e x = it 8 and e y = 2ir/L y . In the charge 
sector, cusps are seen at (2e x ,0) and (0,2e y ). There is 
however another feature namely a cusp at the wave vector 
(e x ,e y ) (Fig. ||b). For this lattice size the dynamical spin, 
S s (q,u>), and charge, S c (q,u>), structure factors are plot- 
ted in Fig. [|c,d. At the locations of peaks in the static 
structure factors soft modes in the dynamics are present. 
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FIG. 5: Ground state spin and charge degrees of freedom 
on a 30 x 12 lattice U/t = 3, S = 0.2. (a), (6): static spin 
and charge structure factors. The bullets denotes q points 
where cusps or peaks are seen. (c),(d): dynamical spin and 
charge structure factor. The dashed line corresponds to twice 
the quasiparticle gap estimated from the density of states of 
Fig. PJc. For clarity we have normalized the data. The weight 
under each curve may be recovered from Fig. (a),(b) since 
&U)S c/s {q,u) =irS c/s (q). 



The lower edge of the spectra along the x-direction is 
plotted in Fig. ||. On this lattice topology, it is possible 
to arrange the holes in two stripes along the x-direction 
(corresponding to an incommensurate filling) or in five 
stripes along the y-direction. If the stripe ordering turns 
out to be long-ranged, one expects one of the two direc- 
tions to be spontaneously chosen in the thermodynamic 
limit. On the other hand if there is no long-range stripe 
order, just short range fluctuations, we expect this state 
to remain stable as the size of the system grows. We note 
that a slave boson approach to the t-J model yields an 
identical pattern of spin and charge modes. For given 
soft spin modes at wave vectors q = (it ± e, n ± e) charge 
modes are generated at wave vectors q — (±2e, ±2e) and 
q = (±e,±e) fll2]| . Our findings support this point of 
view. 

To summarize, we have introduced a new QMC ap- 
proach which allows us to obtain insight into doped Mott 
insulators. It is based on the observation that when 
the number of fermion flavors N — An, a broken spin- 
symmetry Hubbard model may be simulated with no sign 
problem regardless of the lattice topology and band fill- 
ing. More generally, as a function of growing values of 
N the sign problem becomes less and less severe. Our 
approach clearly provides a way to go beyond Gaussian 
fluctuations around a given saddle point. On the other 
hand, it is at present not clear if the extrapolation from 
finite N to N = 2 is justified. We have used this method 
to investigate the doped Mott insulator at N — 4. We 
find a metallic state with no quasiparticles since there is 
a gap or pseudogap in the single particle density of states. 
The low energy excitations are collective spin and charge 
modes. On our largest rectangular topologies they are 
gapless in the spin sector at (tt ± e x ,Tt) and (7r,7r ± e y ) 
and at (±2e x ,0), (0,±2e y ), (0,0) as well (±e x ,±e y ) in 
the charge sector. 
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